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FREE -STREAM DISTURBANCES, CONTINUOUS EIGENFUNCTIONS, 
BOUNDARY-LAYER INSTABILITY AND TRANSITION 


SUMARY 

The research conducted under this project has been directed toward the 
double objectives of providing (1) a rational foundation for the application 
of the linear stability theory of parallel shear flows to transition 
prediction and (2) an explicit method for performing the necessary 

calculations. 

The fundamental discovery upon which our subsequent work is based 
was that the solutions of the linearized, three-dimensional, incompres- 
sible Navier-Stokes equations u,p and the adjoint solutions u,p 
satisfy a "continuity" equation 

|f ♦ V • J - 0 (1) 

« A ^ 

where ' p is a pseudo- energy density (the dot product of u* and u) and 
J is a pseudo -current. This result is derived and discussed in detail 
in Appendix A. 

We next considered (see Appendix B) the expansion of an arbitrary, 
two-dimensional solution of the linearized stream function equation in 
terms of the discrete and continuum eigenfunctions of the Orr-Sommerfeld 
equation in the half-space, y«[0,»): that is, we considered boundary- layer, 

wake, jet or free-shear layer flows. We used equation (1) to derive a 
biorthogonality relation between the solutions of the linearized stream 
function equation and the solutions of the adjoint problem. This is the 
biorthogonality relation for the mixed initial -boundary value problem. 


For the case of temporal stability, we used equation (1) to derive 
the formal solution of the initial value problem as a sum over the discrete 
modes plus an integral over the continuum functions and showed that this 
expansion is complete. We found that the vorticity distribution at the 
initial time is sufficient information to determine the expansion coef- 
ficients and gave explicit formulas to calculate these coefficients. 

For the spatial stability problem, we showed that the continuum has 
four branches. We used equation (1) to derive the spatial biorthogonality 
relation and the formal solution to the boundary value problem. We have 
(sefe Appendix C) also derived the Fourier (in t), Laplace (in x) transform 
solution of the spatial stability problem and used it to show that our 
spatial expansion is complete. 

The boundary conditions for the spatial problem are the Fourier 
transforms, in time, of the stream function and its first three partial 
derivatives with respect to x, evaluated at x = Q. As it stands, this 
formal solution will not give a physically acceptable solution because, 
given an arbitrary variation with y and t at x « 0 of the stream function 
and its first three partial derivatives with respect to x, disturbances 
which lie on all four branches of the continuum will be excited. Therefore, 
as we show in Appendix B, the spatial wave packet will contain, in addition 
to waves propagating toward x * ®, waves propagating upstream from x = « 
and standing waves whose amplitude increases towards x * ®. 

A condition must be imposed that, for x > 0, all propagating dis- 
turbances are traveling in the positive x-direction and all standing waves 
have amplitudes which decay in the positive x-direction. It appears that 
this should be done by requiring that the stream function and its first 
three partial derivatives with respect to x, evaluated at x * 0, be 
orthogonal, using the spatial inner product, to all eigenfunctions on 
branches 2 and 4 of the continuous spectrum. 

It is easy to see that these two orthogonality conditions reduce the 
number of boundary conditions at x * 0 from four to two. This means that, 
for the spatial stability problem, the proper boundary conditions at x = 0 
are the specification of the temporal Fourier transforms of the velocity 
components u and v, for all y. Although these boundary conditions 
were derived from consideration of the continuum eigenfunctions, they 


apply as well to the discrete , Tollmion-Schlicting modes. We have not 
yet carried out a detailed investigation of the implications of imposing 
this orthogonality requirement on the boundary conditions; however, the 
immediate result that the boundary conditions at x • 0 are the specifi- 
cation of the temporal Fourier transforms of u and v for all y appears, 
on physical grounds, to be correct. 

We have presented preliminary numerical results of the application of 
this expansion method at the Fifteenth International Conference on 
Theoretical and Applied Mechanics (Appendix D). We considered the temporal 
stability problem and a simple initial disturbance. We assumed that at 
t ■ 0 the vorticity C was given by 

C ■ C Q e 0 6(y - y Q ) (2) 

a periodic layer of vorticity at a distance y Q from the boundary. 

The stream function is then given by equation (55) of Appendix B, 
and it is easily seen that the expansion coefficients are [from equations 
(56a, b) of Appendix B); 


A n (a) " c o*n Cy o 5 5(0 _ a o> 

(3a) 

A k (a) « C 0 t£Cy 0 ) ^ (a - a Q ) 

(3b) 


The solution of this simple problem, which is in effect the Greens function 
in y of the initial value problem, shows that the amplitudes of discrete, 
Tollmieri'Schlicting modes and the continuum functions are the products of 
the magnitudes of the corresponding adjoint functions, evaluated at y Q , 
the height of the initial disturbance from the boundary and the vortex 
strengths. 

We applied this result to two different flows. The first is a slip 
flow past a bounding plane at y * 0. Although the base flow velocity 
does not vanish at the boundary, we required that the disturbance velocity 
vanish at y * 0. We found (Appendix B) that, because of the simple form 
of the base flow, all the calculations could be carried out analytically 
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and the stream function could be expressed as a finite sum of exponentials 
and error functions. We found that the disturbance retains its identity 
as a periodic array of vortices for all time, but as time increases it 
diffuses, the vortex strength decays, and the centers of the vortices 
drift away from the boundary. 

The second flow we considered is the Blasius boundary layer. The 
velocity s cale wa s taken to be the free-stream speed U Q and the length 
scale was ^vx/ U Q . We chose o ■ 0.179 and R ■ 580.0. At this a and R, 
there are seven discrete Tollmien-Schlicting modes, one of which is 
unstable. We ’numerically calculated the seven eigenfunctions and adjoint 
eigenfunctions and normalized them so that 


< V 


m 


nm 


(4) 


Plots of the amplitude and phase of the normalized eigenfunction and 
adjoint eigenfunction of the seven modes as a function of y, the dimen- 
sionless distance from the boundary, are given in Appendix D. These 
modes are numbered in order of increasing stability with mode 1 the 
unstable mode and mode 7 the most damped mode. 

The amplitude of a mode, say excited by the vortex sheet at 

y * y , is proportional to the amplitude of <J> n evaluated at y Q . It is 
clear from an examination of these figures that when the vortex layer at 
t * 0 is in the inner portion of the boundary layer, say y $ 2.0 (the top 
of the boundary layer is at y » 5.02), there will be a relatively strong 
excitation of the discrete Tollmien-Schlicting waves. Modes 1, 2, and 
3 will have the largest amplitudes, and the higher modes will have sub- 
stantially smaller amplitudes. It is also quite clear that, when the 
initial disturbance is more than about four boundary-layer thicknesses 
from the wall at t a 0, the discrete Tollmien-Schlicting modes excited 
by the disturbance will have extremely small amplitudes. We believe that 
this result is a theoretical explanation of the experimental observation 
of Kachanov, Kozlov, and Levchenko (1978) that vorticity disturbances 
passing above a boundary layer are very inefficient generators of Tollmien- 
Schlicting waves in the boundary layer. 


' CONCLUSIONS 

, We believe that we have created a rational foundation for the appli- 
cation of the linear stability theory of parallel shear flows to transition 
prediction and given an explicit method to carry out the necessary cal- 
culations. We have shown that these expansions are complete. We have 
also carried out some sample calculations which show that a typical 
boundary layer is very sensiti v ? to vorticity disturbances in the inner 
boundary layer, near the critical layer} vorticity disturbances three or 
four boundary-layer thicknesses above the boundary are nearly uncoupled 
from the boundary layer in that the amplitudes of the discrete Tollmien- 
Schlicting waves are an extremely small fraction of the amplitude of 
the disturbance. 

After the completion of this grant we intend to continue these 
calculations. We will continue the calculations of temporal disturbances 
in typical boundary layers and begin calculation of spatial disturbances. 

LITERATURE CITED 

Kachanov, Yer S.} Kozlov, V.V.; and Levchenko, V. Ya, . Occurrence of 
Tollmien-Schlicting Waves in the Boundary Layer Under the Effect 
of External Perturbations. Izvest. Akad. Nauk SSSR; Mekhan. Zhid. 

Gaza, Vol. 5, 1978, p, 85, 



APPENDIX A 


EXPANSIONS IN SPATIAL OR TEMPORAL EIGENMODES OF THE 
LINEARIZED NAVIER- STOKES EQUATION 


Submitted to Journal of Fluid Dynamics 


Expansions in Spatial or Temporal Eigtnmodes of the 
Linearized Navier-Stokes Equation 

by Harold Sal wen 

Department of Physics and Engineering Physics 
Stevens Institute of Technology, Hoboken, N. J. 07030 


The expansion of an arbitrary flow field in terms of the temporal or 
spatial eigenmodes of the linearized Navier-Stokes (LNS) equations for an 
incompressible fluid is developed from a unified perspective. It is £hown that, 
for (v,p) a solution of the LNS equations for a given base flow and (u,q) a 
solution of the corresponding adjoint equations, a scalar "density", 

£(u,v), and a vector "flux", t (u,q,v,p), may be defined such that £ and f 
are bilinear in (u*,q*) and (v,p) and satisfy the "continuity" equation , 
aE'at + v»r * 0. This equation is then used to derive biorthogonality relations 
between the eigenfunctions and adjoint eigenfunctions of the LNS equations for 
a general translationally-invariant problem. In the temporal case, the inner 
product is ffj£dx * ///%>u**v dx which is the natural extension of Schensted's 
Inner product for two-dimensional disturbances and satisfies the requirements 
for an inner product in a Hilbert space. In the spatial case, the "inner 
product" is fff r x dydzdt which is not positive definite. The formal solution 
of the LNS equations is derived, in terms of the eigenfunctions and the initial 
or boundary conditions, for the temporal and spatial cases. It takes the form 
of the evolution of a three- or six-dimensional vector — (v^, v y , v 2 ) in the 
temporal case or (v x , v y , v 2 » 3v^/3X, 3v 2 /3x, p) in the spatial case. 
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V. Introduction 


A few years ago, Grosch and I, after showing that the Orr-Sommerfeld 
equation for unbounded flows such as the Blasius boundary layer possesses both 

I 

temporal and spatial continuous spectra (Grosch and Salwen, 1975, 1978) . set 
out to find the form of the wave-packet expansion for the temporal or spatial 
evolution of the stream function of an arbitrary two-dimensional “infinitesimal" 
disturbance in terms of the corresponding temporal or spatial eigenfunctions. 

We sought to prove a biorthogonality relation between the eigenfunctions 
of the Orr-Sommerfeld equation and of its adjoint and, thereby, to solve for 
the coefficients of the expansion in terms of the inner products of the adjoint 
eigenfunctions with the stream function at the initial time or position. 

This worked out easily in the temporal case, with an inner product equivalent 
to Schensted's (1960) and only minor complications due to the infinite domain 
and continuous spectrum. In the spatial case, on the other hand, we found that 
we didn't know the appropriate inner product and we couldn't find any papers 
dealing with the problem. I therefore undertook the spatial expansion problem 
and, eventually, was rewarded with the result reported here— a unified treat- 
ment of the spatial and temporal expansion problems for solutions of the 
linearized Navier-Stokes (LNS) equations for an incompressible fluid. 

Section 2 is devoted to the derivation of a “continuity" equation which 
is used, in Section 3, in the definition of the inner products and the derivation 
of the biorthogonality relations. These, in turn, are used in Sections 4 and 5 
to derive the formal solutions of the (temporal) initial value problem and 
the (spatial) boundary value problem, respectively. The application of these 
results to two-dimensional disturbances of a boundary layer has been presented 
in a separate paper (Salwen and Grosch, 1980). 
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In order for the formal solutions derived in Sections 4 and 5 to be 
actual solutions of initial and boundary value problems, the eigenfunctions 
used In the expansions must, form complete sets. Not all the eigenfunction 
sets one might want to use have been proven* to be complete but there are, 
by now, proofs of completeness for large classes of temporal eigenfunctions 
for bounded flows (Yudovich, 1965 and DiPrima and Habetler, 1969) and temporal 
(Sal wen and Grosch, 1980) and spatial (Salwen, Kelly, and Grosch, 1980) 
eigenfunctions for unbounded flows. 


2. "Continuity 11 equation 

I start with the LNS equations for an incompressible fluid with a base 
flow 0, 


V* v * 0 


» [ 3T + 9 - f V 1 + U 1 j 

* V A, - |£ • 1 * >.2.3. 

.i* 

and the corresponding adjoint equations , 

■fU 

* <4* _ 

v-u * 0 



3U.. 

w 


- v 


(0*^ ) 


30* 

ax. 





. i 


1,2,3, 


(la) 


(lb) 


(2a) 


(2b) 


+ The complex conjugate, 0*, is used here in order to obtain the correct formal 
expressions. In most applications, 0 will be real, so 0* * 0. 
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For any solutions (v,p) of (1) and (u,1f) of (2), define* 

Edto-j.S*-*. 


Then 




( -P [(tr &,*) • 7 » ( + v, ? -(0 

- p [(V ») • V U, . |f . (I* v ( )' 


+ u 


U|* ? v» - V| ? U|* 


*1 

M * lE_ 4> Ifl* y 


, 1 3X| ■ 3X^ i 


so, with 


f (u.q.v.p) « (| pu*-v) ft + | w ^ [($ U|*) 
- u^*v vj + ijr [u*p + q* v j , 


we get 


|f £(u,v) + 7-f (u,q,v,p) « 0 , 


( 3 ) 


(4) 


(5) 


( 6 ) 


The constant factor, | p , is included in order to emphasize the 
relation between £ and the energy density, j 


1 p v 2 . 


which has the form of a continuity equation relating the time derivative of 
the ''density", € , to the divergence of the "flux", r , For any fixed 
volume V bounded by a surface S , the "continuity" equation (6) may be put 
into integral form, 

3 ^/// £(u.v) dt + U n*r (u,q,v,p) dS ■ 0 . (7) 

V s 

3 * Application to a steady, translationally-invariant base flow . 

Biorthogonalltv relations 

In this section, Equation (7) will be applied to the case in which the 
base flow and boundary conditions are independent of x and t . For all x 
and t, the base flow lJ(y ,^) , disturbance velocity and pressure (v,p) and 
adjoint velocity and pressure (u,q) are assumed to be defined in a closed, 
bounded area, A , of the y»z plane and to satisfy the boundary conditions 

v(x,y,z,t) * 0 , u(x,y,z,t) * 0 for (y,z) e C , (8a, b) 

on the boundary, C, of A. In this case, the temporal and spatial eigenfunctions 
discussed below will form discrete sets. (The extension to an unbounded area 
is not too difficult (see, e.g., Salwen and Grosch, 1980) but it requires the 
relaxation of the boundary condition (8) and the consideration of continuum as 
well as discrete modes.) 

Because of the choice of base flow and boundary conditions, (1) and (2) 
are now invariant with respect to translations in x and t and, therefore, 
possess solutions of the form 

t(x,y,z,t) * t 0 (y,z)e^ oX ' ut ^, (9a) 

U 






p(x»y,z,t) » p 0 (y,z)e^“ x " wt ^, 

(9b) 

fi(x,y,z,t) « u o (y,z)e^ 0x ’ 

(10a) 

q(x,y,z,t) ■ q o (y,z)e^ 0x " 

(10b) 


Because of (8), 3 r 2 * 0 on S. Then evaluation of (7) over a thin slab 

perpendicular to the x-axis for functions of the form (9) and (10) gives 




r x (u,q,v,p) dydz 


» i(e* - a) 11 r x (u,q,v,p) dydz * (11) 

which will be used here to prove biorthogonality relations for the spatial 
and temporal eigenmodes. 


The temporal eigenfunctions are the solutions of (1) and (8a) having the 
form (9) with a real. These may be denoted by (’tf , p ), corresponding to 

c*n otfi 

the x,t variation e^ aX " “n^ 01 )^). For each such solution, there is an adjoint 
eigenfunction (u an q an ) which is a solution of (2) and (8b) having the variation 
g i(ax - v n (a)t) W f th v n (a) = w*(a). Application of (11) to these functions 
gives 

(u m (“) “ “n (o)) |f £ ( ^m * dydz = 0 ’ (12) 

A 

so that the integral vanishes when ui m (ct) f a) n (a) and, with appropriate 
normalization, 


12 


03 ) 


( u i v ) dydz 
'“am* air * 


fur/ 2 ” ' 


(The expression is a function of x and t but is constant because the exponentials 
in the two factors cancel.) This biorthogonality relation for fixed a leads to 
the result for the full set of temporal eigenfunctions, 



V 


- = \ fl V dyd2 dx 

A 


| e i(B ’ a) * dx f| £(S 


\ 

am, v 3iv 


x*Q 


dydz 


2n5(a-tj) || £(S, 


V ) 

am, an' 


X*0 


dydj 


■ «(»-8) ( m for all t . (14) 

The spatial eigenfunctions and adjoint eigenfunctions are the solutions of 
(1), (2), and (8) having the forms (9) and (10) with u and v real. These may 
be denoted by ? un = (v un p wn ), with the variation e 1 ^“ n ^ x " and 

5 vm * (u ym %,„)> with the variation e^ Bm ' v ^ x " As in the temporal case, 

the eigenfunctions and adjoint eigenfunctions may be paired, with B n<“> ■ V<“> 
in this case. The analogous results to (12) and (13) are 


(a (w) - a (u)) r (u o t P ) dydz s 0 

' V 1 n v " jj x vu um, \n, wn,*w * 


and 


<v 


r v (u m q „ v „ p, J dydz = s/2it , 
x ' urn, un , uin/un 7 J mn * 


(15) 


(16) 
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which lead to the biorthogonality relation for the full set of spatial 
eigenfunctions 




( L, 


A* 

‘W* vn , vn 


) dydz dt 



« $(u-v) 3 jJm for all X . (17) 

We thus see that, with appropriate "inner products" <,* and M- 
the temporal and spatial eigenfunctions satisfy biorthogenality relations with 
the related adjoint functions, The temporal inner product <,> satisfies all 
the conditions ordinarily required of an inner product. The spatial inner 
product U , on the other hand, is not positive definite. This is related 
to the fact that disturbances can propagate in both the downstream (+x) and 
upstream (-x) directions. 

4 . Temporal expansion of an arbitrary solution of the INS equations 

The temporal inner product introduced in (14), 

to 

«u* v> * | || u*(x,y»z) • v (x,y,z) dydz dx , (18) 

-® A 
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Is defined for any pair of ordinary vector functions of position. In 
particular, when applied to a solution (v, p) of (1), it involves the 

velocity, v, but not the pressure, p. It Is natural then, in seeking an 

expansion solution in terms of the temporal eigenfunctions, to expand v 
alone in terms of the velocity part {v an > of the eigenfunctions. 

Let (v,p) be a solution o*’ (1), satisfying the boundary conditions (8a). 
Assume that v can be expanded in the form 

CO 03 

v(x,y,z,t) ■ j c on (t) v an (x,y,z,t) da . (19) 

n*l 

Then, by (14), the coefficients are 




**► 

<U 


an , 


v> 


( 20 ) 


so that (using (7) ) 


dc 

an 

dt 


ff ft d * dz dx 


00 

- f \\ ?,f (3 an, P> d » dz d * ■ 0 


and 




CL ( t ) a C „(0) * <u _ ~v> 

v»n' an an, 


t » 0 . 


( 21 ) 


( 22 ) 


IS 


The result is 


00 00 

S(x.y. Z ,t) - 21 f <%,. *>l y.z.t) 

n«l 


dot » 


(23) 


On the assumption that the expansion can be differentiated term-by-term. 


ix| * • p [jr + & * ^ V 1 + V " * ^ u i 1 


oa co 


ZL | ^ v »n, 

n»1 -« t*o ^ 

3V • 

• p L-rr + i! • *\n • ? “il l-do 


flu 00 


£ H. 

n*l -<» 


v> 


t*o 


3p c. 

3X.j dct » 


an 


so that, except for an additive function oft only, 

00 00 

p(x,y,z,t) * ^ | <u,v>| P an (x,y,z,t) da . 
n»l t “° 


(24) 


(25) 


Equations (23) and (25) are the formal solutions for v and p in 
terms of the initial velocity, v(x,y,z,0). 

5. Spatial expansion of an arbitrary solution 

The spatial inner product, J[ ,T[ , introduced in (16), cannot be eval- 
uated in terms of the values of u , q , v , and p at a fixed x because 
r x involves x-derivatives associated with the second derivatives in (lb) and 
(2b). To get around this problem, one can regard the flow field at a given x 
as a 6-vector and make use of the fact that the velocities under consideration 


have vanishing divergence ((la) and (2a)). 


Let e and n be 6-vectors with components 

3V 3V 2 

5 1 “ V h " v y* c 3 * V *4 * V 5 3X^ * *5 * v z' s ax“ * *6 * F 

3U 3U Z 

n l * V n 2 “ u y* n 3 “ V n 4 * u y* 5 ax* - * n 5 * u z' " ax""’ n 6 * q 


and let $ • u * $ • v * 0 . Then, in terms of these components, 


r x (u,q,v,p) » h p(u/v x + u *v + u z *v 2 ) U x 


■x 'x y y 

3V., 3V, 


r 3V u 3V 7 

+ h v^V^ay 1 * aT* ■ v x ( 


3U * 3U * 


ay az 


-) 


+ (v y u y ** - u/v y ') + (v 2 u z ‘* - WQ 
+ (u x *p + q*v x ) 


so that 

r x (u.q.v.p) dydz dt 

may be evaluated in terms of the components of 5 and n at fixed x 
This choice of coordinates also eliminates second derivatives from (1) 

which becomes 



a 


ox 


(30c) 


t 




tills. 

u ox 


v x + 


. 0U_ 

£. — 2- v + 

y oy y 




+ e. 

y 


U v 1 
x 2 


+ ll£ 

U 02 


(30e) 


iL * 
ox 


[ 


0 2 V v 0 2 V ov 

y(au2T“ + g2 2 ” P <0t " + ^ 


3V. 


OV. 


‘Oy* 


Ov 

+ p (u x 5 y 


ou 


y ay 

ov_ 


+ U — — + — — v 1 
u 2 DZ * v ' 


5< 

ox 'x 


j 


oy v y^ + p ^ U x 02 


OU ov ' ov ' 

02“ V 2 ) ‘ P ( l/“ + * (3 ° f) 


It is now straightforward to carry out the formal solution for the 
spatial expansion. The expansion is 


?(x»y, 2 ,t) 


c (J x) 5 un (,t ’ y>2 ' l) *> • 


n«l -« 


with coefficients 


c m M * £ \n, S J- 


(31) 


(32) 
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Then 


dc n 

(jjV 

lx 


m 

| (( 'k r x dydz dt 


3 / •¥ Aj 4 \ 

ay r y ^ u «n*%n, v ‘ p ^ 


id M >4 *4 \ . d 

* 5 ? r 2 (u u n.V' ,) + 5 t 




dydz dt 


,<*> * ‘J<* * (I 


The solution is then 


eU.y.z.t) = 


lv '1 


<=o e jn ()l - y ’ z > t) *■ 
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ABSTRACT 


The expansion of an arbitrary two-dimensional solution of the 
linearized stream function equation in terms of the discrete and continuum 
eigenfunctions of the Orr-Sonnerfeld equation Is discussed for flows In 
the half-space* y c [ 0, •). A recent result of Salwen Is used to derive 
a blorthogonallty relation between the solution of the linearized equation 
for the stream function and the solutions of the adjoint problem. 

For the case of temporal stability, the orthogonality relation 
obtained Is equivalent to that of Schensted for bounded flows. This relation- 
ship is used to carry out the formal solution of the Initial value problem 
for temporal stability. It Is found that the vortlclty of the disturbance at 
t ■ 0 Is the proper initial condition for the temporal stability problem. 
Finally, It Is shown that the set consisting of the discrete eigenmodes 
and continuum eigenfunctions is complete. 

For the spatial stability problem, It Is shown that the continuous 
spectrum of the Orr-Sommerfeld equation contains four branches. The biorthogon- 
allty relation is used to derive the formal solution to the boundary value 
problem of spatial stability. It is shown that the boundary value problem 
of spatial stability requires the stream function and its first three partial 
derivatives with respect to x be specified at x ■ 0 for all t. To be appli- 
cable to practical problems, this solution will require modification to 
eliminate disturbances originating at x * * and travelling upstream to 
x ■ 0. 

For the special case of a constant base flow, the method Is used to 
calculate the evolution in time of a particular Initial disturbance. 
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1. Introduction 


Recent calculations of the discrete eigenmodes of the Qrr-ooamerfeld 
equation (Jordinion, 1971} Mack, 1976; Houston, Cornor, and Boss, 1976; 

a 

Murdoch snd Stewarteon, 1977) have indicated that, for a liven Reynolds 
number snd wave nusber (frequency) , the Orr-Scsserfeld equation for Blssius 
flow has only a finite nunber of discrete temporal (spatial) eigenf unctions. 
Since a finite set of functions cannot be complete, these calculations 
raised the question of hov to expand the stress function of an arbitrary 
disturbance in tens of the noraal nodes. These authors suggested that 
in addition to the finite discrete spectrua, which they found, there is 
a continuous spectrus. 

In Part 1 (Gres oh and Salwen, 1978a), we dealt with the existence of 
the continuous spectrua and the fora of the related eigenfunctions for 
both the temporal and spatial problems . • we showed that the Orr-Sommerf eld 
equation, for any mean shear flow approaching a constant velocity in the 
far field .possesses a continuous spectrum; we gave formulae for the 
location of the temporal and spatial continue in the complex wave-speed 
plane; and we calculated the temporal continuum eigenfunctions for some 
particular cases. In this paper, we turn our attention to the use of the 
discrete and continuum eigenfunctions of the Orr-Sccmerfeld equation to 
calculate the temporal or spatial evolution of an arbitrary solution of 
the linear disturbance equations. 


2a a recent critique of the application of stability theory to the 
prediction of transition, 3ar*ar and Aroaaty (1977) point out that, on 
tha basis of tha liaitad experimental evidence that is available, tha 
coupling of free stream dlsturbancas to disturbancas in tha boundary 
layar appears to be extraordinarily weak and extremely selective in 
frequency and wavenumber, Mack (1977) aakas tha saaa point in a diffarant 
way. Ha points out that "if thara vara no disturbancas Einside tha 
boundary layar], thara would be no transition and tha boundary layar would 
remain laninar. Consequently, it is futila to talk about transition 
without in seme way bringing in tha disturbancas which causa it ... ", 

Mack aids, "... tha pracisa sachanisa by which, say, fraa straaa turbulanca, 
sound, and diffarant typas of roughness causa transition rasains to ba 
discovered. " 

Tha aost datailad discussion of this problaa appaars to ba that of 
Obraaski, Morkovin, and Lindahl (1969). Thay considar various possibla 
mechanisms by which sound or vorticity waves in tha fraa strain night 
intaract with tha boundary layer and cause transition. On tha basis of 
tha available experimental evidence, they conclude that only a snail 
portion of the external disturbance field excites Tollni en-Schlichting 
(T— S ) waves in the boundary layer and a significant portion appears to 
travel within tha boundary layer with little or so interaction. The 
(unstated) conclusion seeas to be that the mechanism which couples free 
stream disturbances to a boundary layer and, thereby, initiates transition 


is unknown. 


Tha cancral problsm bars is tha solution of tha ganaral initial 
and boundai'y value problems for disturbances to boundary layer flow-— how, 
givea the form of the disturbance at a time t ■ 0, to find its variation 
with time and how, given the fora of the disturbance at all tiaes on a 
plane, x ■ 0, perpendicular to the boundary layer, to find out how it 
propagates downstream. In this paper, we approach these problems, in the 
approximation obtained by assuming parallel flow and linearising with 
respect to the disturbances, by expressing the solution as a sum over the 
discrete normal modes plus an integral over the continuum eigenfunctions 
of the Orr-Sommarfald equation. If the (discrete plus continuum) 
eigenfunctions fora a complete sec, this approach will yield a valid 
solution of the problem, 

Starting with Haupc (1912), a number of authors have dealt with the 
completeness of the set of temporal eigenfunctions in a bounded domain. 
Haupc showed that the eigenfunctions for two-dimensional disturbances to 
plane Couecte flow form a complete set and Schensced (1960) proved 
completeness for the eigenfunctions for two-dimensional disturbances to 
plane ?oiseuille flow and for axi-symmetric disturbances to Poiseuille 
flow in a circular pipe. Yudovich (1965) and DiPrima and Habetler (1969) 
have proven completeness of the eigenmodes for a large class of bounded 
flows. We are unaware of any work on the completeness of the spatial 
eigenfunctions or, previous to this paper, on the completeness of the 
temporal eigenfunctions in an unbounded domain. 


In Section 2 we formulate the stability problem for two-dimensional 
disturbances to a parallel shear flow, U(y), 0 < y < », In terms of the 
linearized equation for the stream function and boundary conditions. 

We next formulate the adjoint problem, A new result of Salwen (1979) Is 
then used to derive a pseudocontinuity relation Involving solutions of 
the linearized equation for the stream function and the adjoint solutions. 

This relation Is then used to find the general biorthogonality condition 
for wave-like disturbances to the flow. The biorthogonality relation is 
specialized to the cases of temporal and spatial stability. The orthogonality 
relation for- the temporal stability problem Is that derived by Schensted (1960) 
and discussed by Reid (1965). 

The temporal stability problem is considered in detail in Section 3. 

The solution is Fourier analysed with respect to x. Then the formal solution 
of the initial value problem for the temporal stability of a two-dimensional 
disturbance to a parallel shear flow Is expressed as an expansion in terms 
of the eigenfunctions. The expansion coefficients are determined by inner 
products between the initial disturbance and the eigenfunctions of the adjoint 
equation. We show that the disturbance vortlclty at t * 0 is the proper Initial 
condition for the temporal stability problem. 

In Section 4 we examine the question of the completeness of the set of 
expansion functions for the temporal stability problem. Very recently, 
Gustavsson (1979) has treated the temporal initial value problem by using 
Fourier-Laplace transforms, He finds poles In the transform plane which 
correspond to the discrete T-S modes and a branch cut which corresponds to 


the continuous spectrum. We show In this section that the Fourler-Laplace 
transform solution of Gustavsson is Identical to our Fourier transform, 
eigenfunction expansion solution for the Initial value problem of temporal 
stability. We therefore conclude that our expansion set Is complete. 

The spatial stability problem is considered In detail In Section 5, 

The solution Is Fourier analyzed in t. The formulae for the four branches 
of the continuous spectrum of the spatial stability problem are derived and 
discussed. The formal solution of the boundary value problem for the spatial 
stability of a two-dimensional disturbance to a parallel shear flow Is 
expressed as an expansion in terms of the spatial eigenfunctions. The 
expansion coefficients are determined by Inner products between the boundary 
conditions at x « 0 and the eigenfunctions of the adjoint equation. The 
boundary conditions at x « 0 are discussed. We have not yet been able to 
prove completeness for the set of expansion functions of the spatial stability 
problem. 

In Section 6, we apply the results of Section 3 to the simple case of a 
constant base flow. In this case, we find the eigenfunctions and calculate 
and discuss the temporal evolution of a particular Initial disturbance. 
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2.1 Formulation of the problem 


The besic flow under eoneldereeion is e parallel shear flow, U(y), 
ehe semi-infinite region, y > 0. We are concerned with che temporal 
or spaeial development of an “infinitesimal", two-dimensional distur- 
bance to this flow, (u(x, y , c) , v(x, y, t), 0). In this case, u and 
v can be expressed in terms of a stream function, tP(x, y, t), by 


u • 



( 


v 


'H 


( 


and the linearised Havier-Stokes equations reduce to a single partial 
differencial equation, 

<£ + »£> *** *** •« <■ 





Xn addition, £* must satisfy two boundary conditions at y ■ 0, 



-v(x, 0, t) • 0 


( 5 ) 


and 



u(x, 0, c) • 0, 


( 6 ) 


ad a "f ini tanas s'* condition 



As a consaquanea of aq. (7), iP must satisfy boundary conditions at 
infinity , 






as 


-* + ». 


( 8 ) 


For fixad x and t, ?(x, y, t) balongs to a manifold, M,o£ functions, 
♦(y), satisfying 




2 3 4 

di dfj O ili 

dv* 2* 3’ 4 

ay dy 4 dy J dy 


all dafinad on CO.. ®) , 


(9) 


continuous on [0, ») , 


( 10 ) 



*(0) • 0 and V (0) ■ 0 


(ID 


and 


f* i$(y)| 2 dy 

Jo 


and 


r. »** 


bo eh exise. 


( 12 ) 


The continuum eigenfunctions which will ba discussed in Sections 3 and 5 
do noe satisfy eq. (12). Instead ehay belong eo a manifold M' o M of 
functions satisfying aqs. (9) - (11) and a weakened condition. 


$(y) and bounded in CO, •), 

dy 


(13) 


We define an inner produce, 


(f, g) 5 f f*(y) g(y)dy, (14) 

J o 

in M. The star denotes the complex conjugate. This inner produce is 
defined for the full Hilbert space of functions satisfying eq. (12) and, 
in that space, has the usual properties of inner products. 


2,2 Th« adjoint p rob Ian 


For funceiona £, g c M wa dafina cha adjolne , «C + , of X in eha 
usual wa y by 


{f(x,y,e)}*ff + g(x,y,e)}dxdyde 


■ llj* QT f(x,y,e)}*{g(x,y,e)}dxdyde + Boundary Taras. (15) 


Tha daflnlcion of cha adjolne usad hara yialda an adjolne oparaepr which 
la idaneical eo eha foraal adjolne (Friedman, 1969, pp. 2,3). 

An adjolne aeraaa fuaceion, !? (x, y, c), la a aoluelon of eha 
adjolne aquaclon (vleh U* - U) , 


* +§ ■&+ + + i 


'§ - 0 . 


(16) 


wleh eha boundary condieions at y • 0, 



3f 

ay 


x,0,C 


(17) 


and eha finiteness condition 
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f <!|fl 2 * iffl 2 >“r 


r mi* 

Jo 


+ |*Hdy < • . 


( 18 ) 


Ac above, aquation (19) imp lit* that $ suae satisfy boundary conditions 


a| 3$ 

3x * 3y 


* 0 aa 


y ■* «*. 


(19) 


Whan, as balov, we look for solutions to cha linearised straam 

* 

function aquation (3) which hava a wavalika behavior in x and c, 
aquation (3) raducas to tha Qrr-Somaerfeld aquation and aquation (16) 
raducas to cha adjoint Orr- Sommer raid aquation. Our adjoint Orr- 
Somaerfeld aquation is cha complex conjugata of cha adjoint aquation 
darivad by Schansead (1960) and quotad by P.eid (1965). Tha reason for 
this diffaranca is chat wa dafina cha inner produce in tha usual way, 
(14) while Schanscad's definition of cha inner product (f, g) involves 


f instead of £*. 


2.3 Siorthogonality 


Salwen (1979) haa shown chat eha solutions of tha linearized, 
chraa dimensional Navier-Stokas aquations, u, p and tha adjoint 
solutions u, p satisfy a ''continuity" aquation 


3o . 

1Z + 


0, 


( 20 ) 


whara 


A 

p • u * u , 


( 21 ) 


J - (fr • u)tf + ^ $ COKiJ)^ - uj (V Ui )3 


+ ? p + tip* , 


( 22 ) 


and, as before, cha scar denotes a complex conjugate. 

For Che two dimensional disturbances considered here we will 
introduce two new inner products. Let 3* be any solution of the original 
problem and S? be any solution of the adjoinc problem, chan define 


<$,*> s r 

J 0 


P dy 


r ( i 1*12 + if a?,* 

q 3x 3x 3y 3y ;ay 


(23) 
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*nd 


litfn I f J x iy , 

» 

with J x , the x component of 4, Using (22) and expressing u, p and 3, 
ff la cans of $ and $ , it can ba ahovn that 

Jo R 3x 3 9x 3x 2 3x 2 ** 3x 3 


(24) 


, , 3 2 ff * , q| *3 2 g , 

z 3x3y 3y * z 3y 3x3y J 


-"**£?**#*#* . 
• 2 ff iMir - 


(25) 


Tha fora of chasa ianar products has baan datarainad by the aquations 
for tha stream function and the adjoint stream function. However, we 
can use aquations (23) and (25) to calculate inner products <f, g> and 
If, gH , evaluated at fixed x and t, of any functions f(x, y, t) and 
g(x. y, c). 
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Ie is straightforward co show that <f, g> is dsfinsd for the full 
Hilbert space of functions which satisfy equation (7) and, in chac space, 
has the usual properties of inner produces. On the ocher hand, £ f, gH 
is not positive definite. This is due eo the face chac it is possible eo 
have wavelike solutions to equation (3) which propagate in either ehe upaereaa 
(-x) or downaereaa (+x) direction, 

With ehese definitions it is easy to show thae 

^ <!?,£>+ I ,!?2! "0, (26) 

for any solutions of the original and adjoint problems. 

If £ and $ are wave disturbances of the form 

w 1 * Wo - <»> 

w 1(ax ’ “ e) • <“> 

equation (26) reduces to 

(«• - «*> « 3L • ■ <«' - °*> 1 ■ < 30 > 


This equation may be used co derive biothogonalicy relations for the 
eigenfunctions of both the temporal and spatial stability problems. 


For eh* etaporal stability problea, a is real and given and a* 

•quals a. Xh* orthogonality relation Cor eh* temporal stability problea 
is then 


(a)' - w*) < • 0 . (30) 

so the solutions of eh* eeaporal stability problea and the adjoint 
solutions are orthogonal unless a* • w*, The orthogonality condition, 
equation (30) , can be recognized as being essentially equivalent to 
that derived by Schensted (1960; pg, 27, *q. (2.2.3)), and discuasad 
by Reid (1963). The only difference is that Schensted* s adjoint solution 
is eh* complex conjugate of ours. 

In eh* case of spatial stability, uj is real and given and u' ■ u and 
the orthogonality relation is 

(a* •• a*) I - 0. (31) 

Thus, unless a' ■ a* the spatial eigenfunctions and adjoine eigen- 
function* are orthogonal with ehe inner product defined by equation 
(25). 
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3. Tha Taaporal Stability Froblaa 

3.1 Tha Elganvaluas and Eiganfunctions 

For tha eaoporal stability problaa wo nodify eho finitanass condition, 
aquation (?) , to 


J^^( | Hf| 2 + lf 2, ^,<-. (7') 

This onouros that tho Fouriar intagral expansion of , 

9(x, y, t) - f %(Y, c)o lax da , (32) 

i 1 ^ 

I 

I 

exists. If vt assume chat S^ a of the fora 

l 

r 

^Cy, t) - * a (y)«‘ ia,e , (33) 

r 


i 

i 



than $ la a solution of tha Orr-Somaarfald aquation 


(L* - iaRC(U - c)L - ■ 0 , 

* a dy’ a 


(34) 


with 


and 


c ■ a/a , 


t, .4-a 2 . 


(35) 


(36) 
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.* 

1 - 

\ 

:k 





Similarly us iiium that tha adjoint solution, Sr i alao satiaflaa 

lx 1 

aquation (7 ) thus ansurinf that tha Fouriar lntagral axpanalon of l 

1 

4 

4 

i 

if (x, y, t) » j* <j|(y, t)a iax do (37) | 

if 

if 

axiats. It is assuaad that $ is of tha fora 

i 

£ a (y, t) - , (38) 

with $ tha solution of tha adjoint Orr-Sonaarfald aquation 

ct 

(I l * «*« - '*>■•«, + 2 f * a • 0 . <”) 

vlth 

c* - u*/a. (40) 

Both $ and £ satisfy tha boundary conditions 
a ct 

d a ’(0) - $ 3 (0) 4 3 (o) « i 3 (o) - o , (4i) 



♦a * ♦i * ■” *a * 0 


y 


(42) 


If $ a tad ^ art In M, or cht weaker condition 

♦ a . V ? a , V a faound * d aa y - (43) 

« V 

if $ a and <^ a are in M . Those eigenfunctions which belong to M will be 

called discrete eigenfunctions. Those which belong to M* but not M will 
be called continuum eigenfunctions. 

Xt has been found (Mack, 1976, Grosch and Salwen, 1975, 1978a) chat, 
in general, there is a finite number of discrete eigenfunctions, {^^(y)} 
with eigenvalues {oj^} and a set, ty^}, of continuum eigenfunctions 
with eigenvalues which depend continuously on a real parameter, k, 

in the range CO, •). (Note that the k of this paper is equal to ak of 
Part 1.) 

The number of discrete modes, which we shall denote by N(a), depends 

not only on a but also on R and on the form of U(y) and can, in some cases, 

be zero. The adjoint eigenfunctions also include a finite set, ,. of 

discrete eigenfunctions and a continuum, , with eigenvalues {oj*m> *nd 

{oj*k> » respectively (see discussion following (30)). For a given k, 

0)0 

and vary like a linear combination of e“ 7 as y *. We therefore 
find that 

rk+e , rkK _ , 

j ♦ ct k ,(7)dk j *ak'^ 7 * dk 6 M 

k-e k™€ 


40 


and that, for any square- intagrab la f, 


(f* * ak >. («, ? ak >. <*, %*>. <«. %*> *U txiae. (45) 

Inner produces between continuum functions! such as <J ak » $ k »> do 
noe exist in eha ordinary sansa but ara dafinabla in earns of tha Dirac 
d-f unction (Lighthill, 1960, pp. 10-21), 

Tha discreta eigenvalues aust ba saarchad for (Mach, 1976) , but tha 
continuum eigenvalues follow froa tha asymptotic fora C^ ak » $ ak - linear 
combination of a 4 ^) of eha eigenfunctions as y • » and U * 0^ • U(«): 

(-k 2 - a 2 ) 2 - (ioRU^ - iRw ak )(-k 2 - a 2 ) - 0 , (46a) 


(-k 2 - a 2 ) 2 + (iaRO l - iR«J k )(-k J 


a 2 ) 


(46b) 


so that boeh aquations yield 

u) ok ■ (^) (k 2 + a 2 + ioHUI x ) . (47) 

Ve also find ehat no continuum eigenvalue is also a discrete eigenvalue. 
Than 
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A , > 


0 


for 


u) i* 
cm 


U) 9 

cm 


(48.) 


i 


< ’an' ^ak'* 


%k> ♦, 


an' 


0 , 


(48b) 


«nd 



unlttt ki < k < kj. (48c) 


With proptr labtiling and normalisation, 1c is thtn poasibl* eo chooat 
cha tiganfuncciona in such a way that 



'W 


> • 


V'* 


(49a) 


<A * 0 . ,> * <4 * * $ t> • 0 i 

v an' -ak' v ak* y an' 

and 

4 c*. V” • «< b - *’>• 


(49b) 


(49c) 
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3.2 Expansion of an arbitrary disturbance. 

If tha aiganfunctions form a complete set, then, for any time, t, 
ve may expand *^(y» t) as a linear combination. 




N(a) 

2 

n«l 


*.<“• ‘>* m w * C 


* k (a, t)$ ak (y)dk 


( 50 ) 


of those eigenfunctions. To find the coefficients (a R } and {a^} we may 
make use of eq. (49) to take inner products 


H(a) 

%n' c)> " J ml V (a ’ t)6 nn' * *a (ot ’ 0 ’ 


(51a) 


<$> ak , ^(y, e)> ■ j* a^,(a, e)5(k - k')dk' ■ a^a, t) . 


(51b) 


Ve then find that 


3a (a, t) 

u 

3t 


< ^an* 3t 


3 Ip 
a > 


* iu an %n* V 7 * c)> " ’ iw an V a ’ c) 


(52a) 


And, similarly, 
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5 V°- e) . 

— n *♦<*• ir" ■ - lu ak *k" ,> *>• (!Jb> 

so eh*t 

* (a, c) • A (cOo* 1 ^ 06 (53*) 

n n 

^(a, t) - A k (a)*" iulakt (53b) 

wh*r* 

A ft (a) 5 * n (a, 0) - 4^, ^ a (y, 0)> (34*) 

A k (a) 5 *^(a, 0) - 4 ak , <^(y, 0)> . (54b) 

Th«n, r«£*rring to squacion* (32) and (50), v« find 

y(x, y, e) ■ { S A (a) A (y) «” iu an t 

J-o. a .i « an 

+ A k (a)^ alt (y)*‘ ia:ak,: dk}« iC ‘ x do, (55) 


wh*r« 


A n (0 ° * %n’ ^x (y » 0)> 


■ iff r l r y * °)»* iox ^y 

0 — (36*) 

■ i? £♦»<*> jl C7 ‘- 

and, similarly, 

A k (a> * Iff *ak (y) ] t-0 •’ iaX ^ dy * (56b) 

If eh« discrete and condnuum eigenfunctions fora a conplate sac, 
chan aquation (53) constitutes an expansion of tha scream function of an 
arbitrary disturbance in cents of the discraca (Tollmien-Schlichting) 
and continuum wave solutions, 

, m<7 , , 

of the disturbance equation, (3) , vith coefficients determined by the 
initial form of tha disturbance (x, y, 0). In the next section we will 
show that the discrete and continuum eigenfunctions are a complete set. 
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On* interesting end significant result of this calculation is that th* 
initial distribution of vorticity, 


X 0 <*. •/) - T(x» v> o) - ”|“ - ^ y* °)» (57) 


is sufficient information to determine th* co«ffici*nts (ct) and A^Ca) 
and, therefor*, th* subsequent development of th* disturbance. 


4. Completeness of the Temporal Expansion Functions 

Gustavsson (1979) has carried out a formal solution of the initial value 
problem of temporal stability for three dimensional disturbances. He uses 
the same coordinate system as we do with the addition of the z coordinate 
In the cross stream direction. The formal solutionis obtained by taking 
Fourier transforms in both x and z and a Laplace transform in t , formally 
solving the Orr-Sommerfeld equation in the transform space, and formally 
inverting the transforms. If we eliminate the z-variatlon of Gustavsson's 
solution and his Fourier transform In z (replacing his k by |a|) the two 
solutions should be identical. Both Gustavsson and we express the solution 
in physical space as an inverse Fourier transform over a, the transform 
variable in the x direction. In order to show that these two methods yield 
identical results it is therefore necessary to show that his formal solution 

a ★ 

In Fourier space, v as given in (G13) , is equal to the factor in .curly 
brackets in our equation (55). 

In order to do this we must first translate Gustavsson's notation into 
our notation. Setting g ■ 0, after (G3) It is easily seen that we have the 
following correspondence, 


In order to simplify reference to the equations in Gustavsson's paper 

we will hereafter use the prefix G. Page references are also to Gustavsson's 

paper. 
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This paper 

Gustavsson 


1 v/a 

w 

k 

Id 

Is 

k 

<7 

«i 

1 


In (G3) and thereafter. 

Gustavsson gives the formal solution In Fourier space In equation (G18). 

It consists of a sum of the residue values at the poles plus a contour 
Integral along a branch cut. Using the definitionjof W as the Wronsklan, 
the Dj^glven after (G6)^and the equation (G7), it Is quite straight- 
forward to show that the residue, R , at a pole s y is 

R y = (eV/W) lim C(s-s y ) Ca 1 (s)^(y,s) + a 2 (s)^(y,s)l>. (58) 

Therefore the residue consists of a linear combination of n> 1 and $ 2 
the solutions of the Orr-Sommerfeld equation that approach zero as y + 
i.e. they satisfy (G4) and (42). At s * s y , ^ and $ 2 satisfy the usual 
eigenvalue condition at y « 0 for the discrete modes of the Orr-Sommerfeld 
equation, condition (Gb) (at the bottom of page 1603). This linear combination 
thus satisfies (41). Therefore the residue at s y is proportional to our discrete 
eigenfunction «b (y) with eigenvalue u ,, and 

uv Oiv 

.V • .''“mi* . 
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It Is well known (Coddlngton and Levinson, 1955, p. 101, problem 19) 
that CDj/W]*, the complex conjugates of the functions used in (G6), are 
solutions of the adjoint equation (39). It can be seen from the form of 
(G11) and the definition of our inner product (23) that 

«j ■ <0j\v J ■ 3,4, (59) 

so that the coefficient of $ <y). in the residue is the inner product of 

some solution of the adjoint equation with $(y,Q). Finally, some straight- 
forward, but tedious, algebra shows that the particular linear combination 
of the Oj involved satisfies the boundary conditions (41) and (42) and 

therefore is a multiple of our $ . We thus find that the residue at s is 

av v 

\ ■ ^ V • (60) 

with d av Independent of y ? ( nd A y (a) given by (54a). Before determining d 3tv , 
we turn to the contribution of the branch cut. 

Using the fact that our u * is, it is clear from (G14) that the branch 

cut In the complex s plane is our continuous spectrum in the complex 

w (or c) plane and that the branch point, u ■ 0, corresponds to the limit 

2 

point of our continuous spectrum at c * U-j - i a /R, with U.| ■ 1. The function 
F(a*k; y) in (G18) is, by (G17)and (G19) a linear combination of the soTtrt+ons 
of the Orr-Sommerfeld equation which are, as y ■+• asymptotic to e’ ay , 
e’^y, and e +iky . It can be shown, using (G19), (G20), (G21), and (G22), that 

F(a,k;0) - <|£) • 0 , 

ay y«0 


( 61 ) 


and so F(a,k; y) Is some multiple of our continuum eigenfunction a ak (y). 
Further, It Is obvious that, In (55), 


e 


Kk* 


. f -1«V # -(« 2+ k 2 )t/R 


(62) 


In (G18) with ■ 1. 

Just as for the discrete modes, the Ca*} , v ■ 2,3,4 in (G21) are the 
inner product of some solutions of the adjoint equation with $ 0 . Using the 
definition of the E ffln In (G22) and the definitions of the Ip m « s given 
in the next to last paragraph on page 1604^ some algebra shows that the 
particular linear combination satisfies the boundary conditions at y * 0 
and so the inner product in (G20) is a multiple of the inner product of our 
continuum adjoint, 1> ok (y) , with the initial condition. Therefore, the integral 
term in (G18) is 


I ■ 


d „k A k (o) *ak (y) e ' ( “ akt dl( • 


0 

with A k (a) given by (54b) and d ak independent of y. 
thus takes the form (In our notation) 


(63) 

Gustavsson's result (G18) 


N(a) ilt f 

K (y.t) - 2. d av A y (a) 4 n Jy) e av 


V*1 


r av- 


d ak A k (a) *ak (y) e ’ i “ akt dk • 


(64) 


Both Gustavsson and we may choose our initial condition arbitrarily, 
provided that the various integrals of this function with the adjoint functions 
exist. If we choose the initial condition that <|» a (y,0) is one of 
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the discrete eigenfunctions, say 4 (y), then 

in (55) 

A n^ " *nm 

(65a) 

AjJa) ■ 0 . 

(65b) 

In Gustavsson's formulation (63) we have 

d A (a) ■ 6 
av v v 1 °vm 

(66a) 

A k (a) s 0 . 

(66b) 

We thus see that 

d av 5 1 ' 

(67) 

If we then choose the initial condition 

k+e 

r 

Y 0 (y .0) * J ^ aic . (y) dk' , 

(68) 

a similar argument shows that 

d d< 5 ' • 

(69) 

Substitution of d flv * d ak « 1, (67) and (69), 

makes eq. (64), derived from 


Gustavsson's solution, identical to the curly bracket in our expansion 
solution (55). We have thus shown that the formal solution obtained by 
Gustavsson from the Fourier-Laplace transform is identical, term by term, to 
our formal expansion solution. 


Since any square-1 ntegrable solution possesses a Fourler-laplace 
expansion, we have shown that our expansion (55) Is complete whenever It 
is valid to separate the Fourler-laplace transform solution Into a sum over 
the poles plus an Integral, over the branch cut-* that is, whenever the sum 
over the poles (discrete eigenvalues) converges. This Is, of course, \1so 
the condition for the validity of Gustavsson's solution. 

For the Blaslus boundary layer, the numerical evidence (Mack 1976) 
indicates that, at a given R and a , the number of discrete modes is finite , 
so that the sum over the poles is a finite sum. If this is so, then the 
above condition is certainly satisfied and our expansion functions form a 
complete set. 

We have shown that the Fourler-laplace transform result and the eigen- 
function expansion result are different forms of the same solution of the 
initial value problem to be chosen according to convenience in a particular 
case. The eigenfunction expansion formulation gives explicit formulae 
(54 a,b) to calculate the expansion coefficients. This allows one to calculate 
the amplitudes of the discrete modes (TS modes) and the continuum functions, 
given the initial distribution of vorticity. 


in 


5. Tha spatial stability ptoblaa 

.1 Tha aiganvaluas and aiganfunctiona 

Tha finitsnaas condition, aquation (7), is aodifiad for tha spatial 
stability ptoblaa to 

♦!#*>*•*«-• <n 

This ansuras that tha Fourlar intagral axpansion of tP , 

¥<». y, 0 • f %<*, y)a“ iu,e dw (70) 

axists. If va assuaa that ^ is of tha fora 

y> • v y) * iax ’ (71) 

than <$ is tha solution of tha Orr-Soaaarfald aquation 

u) 

2 

(t? - iRCtaO - w)L - a . • 0 

a a dy z 4> 

with L a givan by (36). 


(72) 


Similarly, va aasusc that tha adjoint solueion, , alao saeisfias 
aquation (7") thus ana u ring that 


$(x, i, e) - I" ? U (X, y)«* l “ e iu 


( 73 ) 


axiats. Wa aaauaa that 


r /„v*ici x 


Vy*. r) ■ Vr>* 


(74) 


Than 3 la tha aolution of the adjoint Orx-Sonmarfald aquation 

(id 


(L^+ IRC (a*U - u)L a * + 2a* • 0 , 


(75) 


Tha boundary condltiona ara 


♦«<o) - V°> 


♦ w (0) - ^(0) • 0 , 


(76) 


and 


*0 aa y -► » , 

Y W Q) T U) 


(77) 


if <? w and ara in M, or 
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V C V K boundtd M y * * • (78) 

if and art in M* , Aa above , cha eigenfunctions which belong to 
M are the discrete eigenfunctions and those that belong to M* but not M 
are the continuum eigenfunctions. 

Jordinson (1971), Corner, Houston, and Ross (1976), and Murdock and 

Stewartson (1977) have shown that there is only a finite set of discrete 

eigenfunctions, ($ r| (y)}, with eigenvalues (ah The set of discrete 

adjoint eigenfunctions, (S }, with eigenvalues (a* } is also finite. 

uin cun 

The number of discrete modes, N(w) depends on R as well as (U and can be sero. 

In part 1 we showed that, in an unbounded domain, the spatial 
stability problem always has a contiguous spectrum. Since then we have 
discovered (Grosch and Salwen, 1978b) , that the spatial continuum of 
Part 1 is only one branch of a four branched spatial continuum. It is 
quite easy to show the exigence of the four branches of the spatial 
continuum. We look for solutions to equations (72) and (75), ^^(y) and 

as MW 

\ik ^> * for m i iv * n which vary lik* •" 7 m y ♦ » (ch« k ustd 

in discussing the spatial continuum in Part 1 is 2/R times the k used here) , 
Noting that, as v * », U * 0^, a constant, and u', u" * 0, we have 

(-a 2 - k 2 )C-a 2 - k 2 - iaRt^ + iwR) - 0 , (79a) 

and 

(-a* 2 - k 2 )(-a* 2 - k 2 + ia*RU 1 -iuR) • 0. (79b) 
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(Note that equations (79a and b) are complex conjugates.) 


Ic is obvious chsc ehsrs iri four roots, {a^}, j ■ 1, ... , 4 with 
snd % 2 6h * roots of 




iRU l°j 


+ k 2 ^ 


iwR 


( 80 ) 


snd 


a 3 • ik , a 4 - -ik (81 a, b) 

Ths eigenvalue a^, th« root of equation (80) with positlvs real part, 

* 

is the continuum eigenvalue discussed in Part 1. As was discussed in 
Part 1, the eigenfunctions of this branch of the spatial continuum 
are waves propagating in the downstream (+x) direction and decaying in 
amplitude as they travel. In the same way it can be shown that a, is 
the eigenvalue of a continuum eigenfunction which is a wave traveling in 
the upstream (-x) direction and decaying as it travels. 

The free scream speed, U^, can be taken to be unity for a boundary 
layer, wake, or free shear flow. In most cases of Interest u/R < < 1. 

It is easy to see that, with ■ 1, and w/R < < 1, 

3 

a l,2 • 2(w/Y) * iRC i a + + (w/R) 2 /Y 3 3 + 0(~) (82) 

R 


with 


frrcj ni*i tw " * 


H 

i 


i 

i 



Y • (1 + 4k 2 /R 2 ) 1/2 . 

(83) j 

;5 

Da fins, as usual, 

Cha phasa spaad Cj by 

j 


Cj ■ ctjui/la^ 2 , 

1 

(84) | 

j 

if 

Than as k ■* 0, 

a x - u + IRC (w 2 + k 2 )/R 2 3, 

f 

i 

(85») j 

■.I 

i 

: f 


a 2 - -w - iRCl + (w 2 + k“)/R"3, 

}i 

(85b) s 


c l : 1 - iC(u 2 + k 2 )/R 2 ]/(w/R), 

(86a) 


c, - - (u>/R) 2 + i(ui/R)Cl + (u 2 + k 2 )/R 2 3. 

* 

(86b) 

Whila as k ■* », 

« l ” (MR/ 2k) + ik , 

■f 

:] 

(87a) j 

| 


a 2 : ’ a l * 

(87b) 1 


c x - w 2 R/2k 3 - iu/k , 

(88a) f 

: 1 

i 

a 

C 2 : " C 1 • 

(88b) ! 

i 
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Tha damping rate, for eha spatial eigenfunctions, is Im(cs) and the 
phasa speed is &a(c). Tha aquations given abova show that eha eigsnfunceionp 
on branch two of eha spaeial continuum, for boundary layars, wakas, and 
fra* shaar flows, always -hawa boeh a vary large damping raea and a vary 
small phasa spaad. This is in oarkad concrasc to ehosa of branch one, 
which, as was shown in Parc l, or can ba aaan from eha abova rasules, 
contains lighcly dampad aiganfuneciona soma of which hava a vary slow 
phasa spaad and soma of which hava a phasa spaad naarly aqual to eha free 
stream speed. 

Tha spaeial continuum aiganf unctions of branchas 3 and 4 art standing 
wavas in x because ehay vary like 


.ia 3 x . 


# ia 4 x 


+kx 


( 89a ,b ) 


As in the temporal case, eha inner produces between the spatial 
continuum eigenfunctions do not exist in tha ordinary sanse but can ba 
defined as 5 functions, Then, with proper labeling and normalisation, 
it is possible co choose Che eigenfunctions such thac (with Che super- 
script i or j indicating the branch of tha continuum) 


(T A , <4 1 3 ■ C i , 

Y um Y tun nn 

(90a) 

■ 1 ^ ■ W 1 • 

(90b) 


(90c) 
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(T? , i .32 (Tc ' 1( V 

y um* un J n R »*» 


do* dO » 

+ ° u „ 2 - V + 2 -if if- 3 


+ uC (a* » + a i a + a 2)0* i ' 

' wn wn uin wn^ tun ^wn 


d5* dd » d 2 0* 

, , JW y an + um k , }}dy , 

* *■ dy dy dy 'urn 7 


and char* ara analogous expressions for cha inner produces in (90b) 
and (90c) . 










^(w, x) 

— X — 




1 W IT 3 " io ik )(I ik’ 
• ia^ «j[ 4> («, x) , 


(94b) 


so that 


a (w, x) - A (w)a i0 ^ :< , 

a a 


(95a) 


l£ 1 ) (u», x) - A,f i) (u)a ia ^ )x , 


(95b) 


whtrs 


A q (u) X a a (w, 0) - (T? wn , W u (0, y)3 , (96a) 

^(u) Sa^Cw, 0) - I $$ , %(0, y)Z3 . (96b) 

From equations (73), (92) • and (95), we have the formal solution to 
the spatial stability problem for the two dimensional, linearised 
Navier-Stokes equations 

_ M(w) . „ 

r<x. y, t) - f { E A(u»4 n (y)e ia ^ JC 
n- l a ^ 

+ Z f A^ 1) (u)®^ ) (y)e ia ^ >x dk}e”^ ut dw . (97) 
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* 


and chart la a similar expression for (u) . 

This is cha formal solution of eha spatial stability problem 
for an arbitrarily imposed boundary condition at x ■ 0. Tha boundary 
conditions which oust ba spaclfiad art cha Fourlar transforms in tlaa, 
of tha strata function and its first thraa partial derivatives with 
raspact to x, avaluatad at x ■ 0. 

As it stands, this formal solution will not give a physically 
accaptabla solution because, given an arbitrary ?P(0, y, t) and darivativas, 
disturbaneas which lia on all four branchas of tha continuum will ba 
excited. Therefore the solution will contain, in addition to cha vavas 
propagating towards x ■ * and cha standing waves whose amplitude decays 
towards x ■ ®, waves propagating upstream from x ■ * and standing waves 
whose amplitude increases towards x ■ ». 

A condition must ba imposed that, for x > 0, all propagating dis- 
turbances are traveling in tha positive x direction and all standing 
waves have amplitudes which decay in Che positive x direction. It appears that this 
Should be done by requiring that 'if (0, y, t) and its first three partial 
derivatives with respect to x be orthogonal to all eigenfunctions on 
branches 2 and 4 of cha continuous spectrum but we have not yet investigated 
the implications of imposing this condition on tha disturbance stream 
function at x • 0. 
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6. Application co ch« Temporal Development of a Modal Flow 


Za this aaceion, w« apply cha rasules of saeeion 3. co eha simple 
bast flow* 


U(y) • Uj_ • conacanc, y > 0 , (TOO) 

which is a slip flew past a bounding plana at y ■ 0. Though eha base 
flow valociey doas not vanish ae eha boundary* wa still raquira eha 
diseurbanca valociey to ba saro ac y ■ 0. Because of eha simplicity of 
eha basa flow, eha expansion functions ara alaaancary functions. 

In 6.1, we find the expansion functions. In this case, there are no 
discrete eigenmodes; all of the eigenfunctions are continuum functions. 

In Suosaceion 6.2, wa solva for eha elaa-duvalopmane of a pareicular 
inielal diseurbanca by expanding in earns of ehesa eigenfunctions . The 
inicial disturbance chosen is a periodic layer of vorticiey confined eo a 
plana parallel to eha (y ■ 0) boundary. Because of eha simple form of the 
initial diseurbanca and eha simplicity of eha basa flow, ie is possible eo 
obeain eha solution in closed form in earns of error funceions. 

6.1 The Eigenfunctions 

For eha base flow of Squacion(lOO) eha diffaraneial aquation, (34), 
for eha expansion functions becomes 

2 2 

(-^T - a 2 - iotR (U, - c))(-^ - a})i • 0 , (101) 

dy* 1 dy 4 
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with eh* general solueion (for a i 0, c i U^), 


$ • A* 


T1«I y 


+ Be 


i«iy 


+ C.” + De 


-py 


( 102 ) 


vh*r* 


p 2 - a 2 + iaRO^ - c) , (103) 

(In chi* cast of constant U, d must saeisfy ch# Sana diffarancial aquation) . 
In addition, i oust sacisfy Sqs, (11) and (13); i.a., d and d' oust 
vanish at y ■ 9 and be bounded in [0, •). Sine* J a ^ y is unboundad, 

B • 0. To sacisfy eh* boundary condition ac ch* origin, w* ausc eh«n hav* 

d ■ AC*- |o|jr - cosh py + slnh py] , (104) 

P 

which is unbounded as y * • unless p is purely imaginary. Th« solutions 
arc than given by 


p ■ ik; 0 < k < - , (105) 

« ak ■ -i(a 2 + k 2 + iaau I )/R , (47) 

d ak (y) * *ak^ " A ak C(«"' a ^ y - cos l<y) +4^1 sin ky], (106) 
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vhsrs tht norsslizscion constinc , 


A «k * \J1 ’ < ,07) 

Is dsescminsd by chs condition 

<£ ak » $ ak »> " $(k - k'). (108) 

In this case, where the and $ ak are known explicitly, one may 
show directly that, for F(y) any continuous, differentiable, square- 
Integrable function In CO, 58 ), 

m 

' ^ck* F> ^k (y)dk " F(y) ‘ e ‘ ia|y F (°> 009) 

o 

thus confirming that the set of Is complete for functions In M, 

with F(0) « 0. 
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6.2 Th« Taaporal Evolution of An Initial Diaturbanca 

In order to deaonstrate the application of this expansion technique, 
we consider the particular initial disturbance 

i 

-7 2 !P(x, 7, 0) - T(x, y, 0) ■ T 0 « 1:ao * Hy - y Q )* (110) 

a periodic layer of vorticlty at a distance yg from the boundary, 
lollowing section 3.2, we find that the stress function at any tine will 
be given by 


*?(x, y, t) 


r r 


4 0 


L Y a >W*> 


aa^t 


dk e‘ 


da 


on) 


where 


A^(a) ■ TF J Q ? alt ( y> J" T(x, y, 0)e“ iax dxdy. (H2) 

It is easily seen, by substituting Sq. (110) into Eq. (112), that 


V a) ’ W y 0 } 5(ct ' a 0 )( 


(113) 


so that 


£(*, y. e) - jT W < a - a 0 )^ k (y) e" 1 "'* 6 dk e l0x da 

- s 1 * 0 ^ “ Ult) •* < ^ C/R ^(y 0 ) ^(y> •’ k2t/R dk. 014) 


After using Eqs. (106) and (107 ) , for <p^, we find chac each cena in the 
integral is expressible as suds of arror functions. Tha rasulcs art givan 
in an Appendix. From thsss results, it can be shown that, for t + ® with 
y fixed, 


9 



4 -c^ 2 t/R 


cos cy[x - U^t) • (function of y) 


(115) 


and, for y • with t fixed, 

% - a”^ y + cos Cjj(x - U^tt) • (function of c) . (116) 


It is clear chac, even though the individual eigenfunctions used in the 
expansion oscillate with constant amplitude as y •**, the wave packet behaves 
like •“V as y * *. 
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Figura 1 shows contour plots of th« scrsam function for tht disturbance, 

In a frame of raftranca moving with the frea straam velocity, at sin diffarant 
timas. Ma hava chosen ■ 1,0 and y Q ■ 1,0 for tha example shown hart. 
Coneours of tha dlsturbanca straam function hava also baan calculated 
for ocher combinations of values of a Q and y Q and, for these ocher values, 
tha evolution of tha disturbance in time Is quite similar to chat shown 
in figure 1. 

In figure 1 the (+) and (-) indicate the position of the maximum and 
minimum values of the stream function, these maximum and minimum values 
are given in the caption to the figure. Tha flow is counter-clockwise 
around a maximum (+) and clockwise around a minimum (-) . 

It is clear from this figure that the discurbancs, which is a periodic 
vortax sheet at t ■ 0, retains its Identity as a periodic array for all 
time, but as clae Increases it diffuses, the strength decays, and the 
centers of the vortices drift away from the boundary at y * 0. 

We could, of course, generalize this modtl problem by considering an 
initial vorticty distribution In tha y direction. We have not carried out 
this calculation becausa cur intant in solving this modal problem was to 
illustrate the expansion procedure and wa do not chink that it warrants 
furthar elaboration. 
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Figure Caption 

Figure 1: Contours of the disturbance stream function for the model 


problem in a frame of reference moving with the free stream 

velocity at six different times. In this exampley^ ■ 1.0, 

<Xq » 1.0, and yg • 1.0. There are twenty contour lines on 

each plot. The values of on these contours are 0.95 ^max, 

0.85 P max, ... , -0.95 P max. The (+) and (-) indicate 

Che positions where P • P. and . . Note that 

& . - - W 

x min max 


(a) 

(b) 

(c) 

(d) 

(e) 

(f) 


t/R - 10" J 
t/R - 10" 2 
t/R - 10" 1 
t/R - 1.0 
t/R - 5.0 
t/R - 10.0 


. « 


max 


- % 

. v. 

. <p 
• % 


max 


max 


max 


max 


max 


- 0.425 

- 0.359 

- 0.212 

* 0.223 x 10' 
- 0.108 x 10 
■ 0.383 x 10 
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Appendix: Solution of ch« Model Problem 

In section 6 we showed Chat Che stream function for the model 
problem is, equation ( 114 ) > 

<f(x, y, t) • e ia ° (x " UlC) e“ a 0 c/R e' k c/R ^ k (y 0 )^(y)<ik, 


where $ . (y) is given by ( 106 ) and ( 107 ). Substituting for $ . (y) and 

c b Q 

(y Q ) in this integral it is straightforward to show chat, with 
0 

T ■ t/R, 


iPCx, y, t) • e io ° <5t ‘ UlC) {l 1 (a 0 , t) - e a 0 y I 2 (a Q , t, y 0 ) 


♦ a ft «"°0 y I 3 (a 0 , t, y Q ) - e _a ° y 0 I,(a n , r, y) 


2 0 * 


+ 2 W T * y + yJ + ? T * y - 


0' - 2 fc 2 O’ 


- a Q i 3 (ct 0 , t, y + y n ) + a n e’ a ° y 0 l,(a Ql r, y) 


O' 0 


*3 V ^Q 1 



+ t aj i,(a n , t, y - y 0 ) - \ * 2 0 VV f. y + yn )} > 


(Al) 
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" ' '■ ***** ******* i fifin 


where Che functions I. art given by 

j 


t.Co.T) . f X 1 f •’ k2 ' c . k * ]<k 

1 * *0 <k 2 +a 2 ) 2 


cf l C<4 + a 2 t) erfc (at* 4 ) - ^ e^S , 

imT 


(A2) 


I 2 (a, r, Z) i | e' a 


T f* ? 3 cos kZdk 

J .0 (k z +a z ) 2 


• C(1 + 2a 2 t - aZ)e“ aZ eric (at* 4 - aZ/Zat* 4 ) 


+ (1 + 2a 2 t + aZ)e aZ erfc (at* 4 + aZ/Zax*) 


Aar* 4 -a 2 t -a 2 Z 2 /4a 2 t 1 


(A3) 


t 3 (a, T, Z) i I e- a T 


e" k T [ ~ y k -' T - 7 ] sin kZdk 
0 (k + a ) 


— ^ C(2a 2 t + aZ)e aZ erfc (at* 4 + aZ/Zat* 4 ) 
4a 


- (2a 2 t - aZJe^ 2 erfc (at* 4 - aZ/Zat* 4 )} , 


(A4) 
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1 


■] sos kZdk 


I 4 («, T, Z) 2 | t“ a2t [" s'** 1 0 


(k 2 + a 2 ) 2 


■ ~ C(1 - 2a 2 T + aZ)s‘ aZ srfc (at* - aZ/ZaT 5 ) 
4a J 


+ (l - 2a 2 t - oZ)* 02 *rfc (at 1 * + aZ/2ax\ , 


+ 4ai^ # -a 2 t r a 2 Z 2 /4a 2 T ] 

JS 


(A3) 


and, as usual 


erfc(Z) 


_ 2 _ 

/? 



(A6) 
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over the discrete spatial n,odes and branch cut contributions equivalent to 
integrals over the sp atial continua, The two methods are therefore equivalent 
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ABSTRACT 

In this paper we give the solution of the boundary- layer receptivity problem: 
that of determining the amplitudes of the Tollmien-Schlichting modes and 
continuum eigenfunctions of a boundary layer given the form of the velocity 
profile and the disturbance, within the context of incompressible, linear 
stability theory for a parallel shear flow. We give the formal solution 
to the initial value problem for temporal stability and give the proper 
initial condition for this problem. The formal solution of the spatial 
stability problem is also given and the proper boundary conditions at x * 0 
and radiation conditions at x ■ « are discussed. We give examples of the 
application of this method to the calculation of the temporal evolution of 
a particular disturbance in two flows, a constant base flow and the Blasius 
boundary layer. 


*This work was supported, in part, by the National Aeronautics and Space 
Administration under Grants NSG 1618 and 1619. 
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SUMMARY 


The last ten years has seen an increasing use of the theory of hydro- 
dynamic stability to predict transition in boundary layers. Mack (1977) 
gives an excellent, up to date review of various transition prediction * 

methods. All of these methods include at least one unknown parameter 
A q, the initial amplitude of the disturbance in the boundary layer. 

There are numerous discussions of the boundary- layer receptivity problem, 
that is, the problem of determining Ao given the velocity profile of the 
boundary layer and the disturbance (Obreinski. Morkovin, and Landahl, 1969; 
Mack, 1977; Berger and Aroesty, 1977). All of these authors conclude 
that the mechanism by which free-stream vorticity and sound disturbances 
generate Tollmien-Schlichting waves in a boundary layer is unknown. 

In this paper we give the solution of the boundary- layer recep- 
tivity problem within the context of incompressible, linear stability 
theory for a parallel shear flow. The expansion of an arbitrary two- 
dimensional solution of the linearized stream function equation in terms 
of the discrete and continuum eigenfunctions of the Qrr-Sommerfeld 
equation is discussed for flows in the half-space, y e[0, ®) . A recent 
result of Salwen is used to derive a biorthogonality relation between the 
solution of the linearized equation for the stream function and the 
solution of the adjoint problem. 

For the case of temporal stability, the orthogonality relation 
obtained is equivalent to that of Schensted (1960) for bounded flows. 

This relationship is used to carry out the formal solution of the 
initial value problem for temporal stability. It is shown that the 
vorticity of the disturbance at t = 0 is the proper initial condition 
for the temporal stability problem. 

For the spatial stability problem it is shown that the continuous 
spectrum of the Orr-Sommerfeld equation contains four branches. The 
modes on these brances are (1) waves propagating downstream, (2) waves 
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propagating upstream, (3) standing waves whose amplitudes decrease down- 
stream, and (4) standing waves whose amplitudes decrease in the upstream 
direction. The biorthogonality relation is used to derive the formal 
solution to the boundary value problem of spatial stability. We show 
that the boundary value problem of spatial stability requires the stream 
function and its first three partial derivatives with respect to <x be 
specified at x * 0 for all time. The imposition of a radiation condition 
downstream* i.e. at x * «, eliminates disturbances which originate at 
X * ® and travel upstream to x * 0, The imposition of this radiation 
condition reduces the number of independent boundary conditions at x * 0 
from four to two. 

We give two examples of the application of this method to calculate 
the temporal receptivity of boundary layers to a disturbance. We specify 
the disturbance at t * 0 to be a vortex sheet parallel to the boundary and 
sinusoidal in the streamwise direction. We then calculate the evolution 
in time of this disturbance in Cl) a constant base flow, for which the 
calculation can be carried out analytically and (2) in the Blasius boundary 
layer for which we calculate the amplitudes of the discrete Tollmien- 
Schlichting waves and of the continuum eigenfunctions numerically. 
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